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Abstract 



This paper presents a methodology and software tools for parametric design of 
complex architectural objects, called digital or algorithmic forms. In order to 
provide a flexible tool, the proposed design philosophy involves two open source 

^ utilities Donkey and MIDAS written in Grasshopper algorithm editor and C++, 

respectively, that are to be linked with a scripting-based architectural modellers 

P y Rhinoceros, IntelliCAD and the open source Finite Element solver OOFEM. The 

r \ emphasis is put on the mechanical response in order to provide architects with a 

• consistent learning framework and an insight into structural behaviour of designed 

O objects. As demonstrated on three case studies, the proposed modular solution is 

capable of handling objects of considerable structural complexity, thereby acceler- 

^ ating the process of flnding procedural design parameters from orders of weeks to 

Q\ days or hours. 

o> 

'^ Keywords: Algorithmic design; Procedural design parameters; Conceptual 

^ phase; FE analysis; monolithic versus modular solution 

cn 
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^_^ 1. Introduction 

• • 

.J^ Architecture, the essential component to urban environment, impacts every in- 

/\ dividual, which, in turn, influences back the next generation of art. In each period 

^ of civilisation, architecture has reflected the level of societal progress by integrating 

the state of the art from various flelds of human activity. In other words, we can 
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a) b) 

Figure 1: Comparison of a) architrave and b) amorphous forms. 

understand architecture as a multidisciplinary subject combining current knowl- 
edge not only from technical fields but also from Humanities, Ecology or Military 
and defence. However, the increasing level of knowledge characterised by narrow 
specialisation results in educational institutions producing architects unprepared 
for a strong cross-disciplinary dialogue vital in today's complex society [Ij. 

The lack of discussion and mutual understanding is evident especially between 
architects and structural engineers as mentioned in the work of Clune et al. j2]. 
As long as the architects traditional, especially architrave, structures where the di- 
mensions of particular components were the only unknowns, see Fig.[T^, a complete 
structural assessment could be performed at late stages of the design process. Cur- 
rent designers, however, often employ sophisticated computer aided environments 
to generate complex amorphous light-weight forms, thereby requiring a conceptual 
static assessment already at the beginning of the design process, see Fig. ^. 

1.1. BIM concept 

To start our discussion on interdisciplinary cooperation, let us first recall the 
Building Information Modelling (BIM) concept [3j. BIM is a recent and popular 
way of managing complex collaboration and communication processes among ar- 
chitects, structural engineers and construction industry members. The term BIM 
involves the process of generating and managing building data throughout the 
life cycle of a structure. The result is a data-rich, object-based, usually three- 
dimensional "Building Information Model" created by specialised CAD-BIM sys- 
tems. It integrates all the information on the construction from architectural de- 
sign (geometry of building elements, spatial relations as connectivity, etc.), struc- 
tural design (project design documentation, structural scheme) to the process of 
construction and maintenance (detailed design, building process and/or rehabili- 
tation). Thanks to this, architects and structural engineers (and also builders and 
owners) can effectively generate and coordinate complex digital documentation of 
the structure at any phase of its existence. 
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Figure 2: Illustration of a) BIM and b) multidisciplinary approach (courtesy of Tuba Kocaturk: 
BIM conference in Prague). 

Despite the obvious advantages, BIM only connects participants of the building 
industry by means of a database-hke communication channel, Fig. [2^. Each par- 
ticipant, however, remains highly specialised in his own field. This is inconsistent 
with our aim to enhance multidisciplinary approach in the design process, where 
the integration of professions into architecture should follow from exchange of mu- 
tual knowledge as sketched in Fig. [2)3. Moreover the BIM approach is in certain 
cases too cumbersome. This happens namely in initial stages of the project (in- 
vestor's plan, architectural study), which is the most creative phase of the design 
process, taking place in a close cooperation between an architect and his client. In 
this case, BIM is unnecessarily complicated and general. On the other hand, this 
phase can last for several months (even years in extreme) and involve considerable 
costs. It is therefore desirable to validate the starting form at minimal time, while 
avoiding severe violations of structural principles. 



2. Methodology 

The methodology and software tools presented in this paper aim to improve 
the interaction between designers and structural engineers in the critical phase of 
the conceptual design. At this time, architects are sorting out preliminary visions 
penetrated by investor's plans. Resulting functional and spatial contexts may be 
difficult to understand even for the members of the comunity, structural engineers 
let alone. In these regards, we identify three goals summarized next. 
(i) collaboration: The developed interface can be understood as a generic tool 
which combines geometric modellers and a software for structural analysis [4j. A 
significant emphasis is given to the modular approach that enables the connection 
among arbitrary open source and commercial software packages. This strategy 
significantly broadens the applicability of each single module, namely, in compari- 



son with recently developed products based on a monolithic solution, e.g. [5l[6l[7j. 
In addition, the set of our tools is released under public license regulations and is 
freely availablq^to corporate and non-profit bodies. 

(ii) learn: From the viewpoint of a designer, the tools are integrated into his 
favourite modeller as a plug-in. It allows for structural analyses of different diffi- 
culty up to the complexity of a target artwork. Probably most importantly, the 
basic interface (GUI) has to be easy to use in order to not discourage a user at the 
first impression. As a result, the software allows the user to understand what he 
does rather then to provide him with plain answers on static admissibility of the 
structure. 

(iii) form-finding: In the case of computationally less demanding structures, 
the analysis runs interactively. The response of the model to loading or geome- 
try changes is visualised in the real time. This, in combination with procedural 
modelling, enables relatively fast generation of a large number of variants and 
instant structural assessment for intuitive shaping of the structure. If necessary, 
such process can be automated by Evolutionary Structural Optimisation (ESO) 
optimisation methods, see ^i9j. 

2.1. Object- oriented model 

As indicated above, there is a fundamental incompatibility in cooperation be- 
tween the architects and structural engineers in terms of priorities imposed on the 
computer model of designed objects. While architects emphasise the aesthetics 
aspects, structural engineers give the focus on the load-carrying system. An anal- 
ysis directly performed on architectural models seams to be the most direct way. 
However, the complex three-dimensional CAD data are often computationally pro- 
hibitive. Moreover, a comprehensive analysis on somehow provisional inputs may 
easily come out uneconomic. It is thus desirable to simplify these models, while 
maintaining their essential structural characteristics. Typically, such a transfor- 
mation is performed by a structural engineer on the basis of his experience and 
professional knowledge, since the full automation of the process is very difficult 
even with the conversion techniques developed within BIM technology. 

As mentioned above our focus is on modelling of preliminary layouts. Thus, 
architects should tolerate simplified models, representing only the load-bearing 
components of the structure. If so, the conversion into the computational model 
can be carried out directly in the geometric modeller with only a minimum expert 
intervention. In the parametric modellers, this is best achieved by exploiting their 
inherent scripting capabilities. 

The conversion can be briefiy illustrated by the example of a straight beam 



^ www. igend.cz 




self-weight 





b) c) 

Figure 3: Object modelling of beam with rectangular cross-section, a) architectural model, b) 
structural model, c) object model. 

with rectangular cross-section, Fig. |3| In the usual architectural model, a beam is 
displayed on the output device and maintained in computer memory as a set of 
twelve lines topologically linked with the nodes located in eight vertices. Fig. |3^. 
For the purpose of an effective structural analysis, this model is simplified to an 
one-dimensional line segment. Afterwards, the computational model is supplied 
with additional information, here e.g. cross-sectional characteristics, material pa- 
rameters and applied loads, Fig.[3]3. A similar object-based approach is also applied 
for planar and shell entities. 



3. Software architecture 

The basic structure of the proposed interface is briefiy outlined in Fig. [4j As 
mentioned above, we exploit a modular approach in which each of the module 
is responsible for a particular action within the communication chain between 
structural engineers and designers. The converter and the plug-ins to geometric 
modellers were newly created (dashed line grey boxes in the fiowchart) and released 
under the open source licence regulations. Existing components were used and 
extended when needed. Where possible, free and open source variants of particular 
modules (round corner boxes) were preferred for their fiexibility and accessibility. 

A Multifunctional Interface Between Design and Mechanical Response Solver 
(MIDAS) [lO] is in the heart of the reported system. It is responsible for manipulat- 
ing input and output data of structural analysis in various formats and was tested 
in combination with in-house developed packages OOFEM [HI |T2], SIFEL ^13] 
and proprietary system ANSYS. As for the input data, there exist several ways 
of generating structural models. For instance, simple benchmarks can be written 
directly in a text editor. On the contrary, unique models are best to be prepared 
by single-purpose generators, see Section [5^3} In most cases, however, the designer 
is expected to come in a close contact only with his favourite modeller and the 
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Figure 4: Flowchart of individual program components. 

corresponding plug-in, e.g. DONKEY fT^^ T!5] and STRUCT ^16j. The remaining 
process is assumed as an automated black-box tool. The plug-in should help the 
user to create a structural model and provide it with additional information to run 
the analysis, see Section [2TJ The sequence of individual routines is as follows: 

1. user: architectural/geometric modeller - generation of structure's geometry; 

2. plug-in: completion of object model; forward VTK export; 

3. MIDAS: data modification and consistency check; generation of finite ele- 
ment (FE) package input file; 

4. FE package: structural analysis; 

5. MIDAS: output data processing; backward VTK export for visualisation 
purposes; 

6. plug-in: visualisation of results. 

As discussed above, typical user is expected to have only a minor experience 
with theoretical and computational aspects of structural analysis and the software 
tool should provide him with an interactive learning interface. To guarantee this, 
the plug-in (s) should meet the following criteria: 

1. only basic structural analysis (linear statics with truss, beam and shell ele- 
ments) ; 

2. geometric model is clean, no confusing details are contained; 

3. material characteristics and boundary conditions can be set up in a simplified 
and extended regime (e.g. predefined or custom materials); 

4. interpretation of mechanical response with optional level of detail that en- 
ables designers to choose a post-processing mode adequate to their particular 
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needs and knowledge (e.g. "yes-no" binary markers, cross-section resistance 
ratio or a full set of internal forces, displacements, strains and stresses); 
5. interactive and intuitive handling. 

The flow of data proceeds through the utility chain by means of flies in various 



formats, see Section |4.l[ since particular modules support different input/output. 
The ASCII VTK (Visualization Tool Kit )q] has been chosen as the primary format. 
It has a human readable syntax and can be visualised directly in the modeller or 
free visualisation tool-kits such as Paraview or MayaVJj Thanks to this, the 
data exchange can be simply controlled at any stage of software development and 
debugging. 

Regarding the FE discretization of geometric models, we have explored two 
equivalent methods. Namely, the modeller triangulation toolkit, originally involved 
for rendering visualisation purposes, and an external mesh generator that is called 
from MIDAS. Since the modern architectural models mostly consist of NURBS 
(Non-Uniform Rational B-Spline) entities native for Rhinoceros [17j, the same 
geometry deflnitions are also used for mesh generation. However, this together 
with built-in generator sometimes leads to poor mesh quality. Thus, a more flexible 
way appears to consist from passing the solid geometry to MIDAS and generate 
the mesh by an external utility, e.g. T3D fT8] or Gmsh fT9l [20j. 

4. Prototype implementation 

The eflicient basis of a implementation is composed of MIDAS and other two in- 
house developed software packages OOFEM and T3D. OOFEM is a modular flnite 
element code for solving problems of solid, transport and fluid mechanics. T3D is a 
mesh generator operating on complex two- and three-dimensional domains. Both 
OOFEM and T3D are compiled in a minimum required conflguration as dynamic 
libraries and linked with MIDAS. The result is released as the open source software 
operating on various platforms. 

4.1. MIDAS 

The module MIDAS [lO] is a tool without graphical user interface designated 
for manipulating both input and output data of structural analysis. MIDAS's 
source code, written in C++, is released under GPLv3+|^ license regulations. It 
can work with data flies of different formats - UNV, VTK, VTK XML as well as 
OOFEM, SIFEL, T3D and ANSYS native formats. 



^ www. vtk.org/VTK/img/file-formats.pdf 

^ www .paraview . org , mayavi . sour cefor ge . net 
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^GPLv3+: GNU GPL version 3 or later, http: //gnu. org/licenses/gpl. html 



Recall that the input geometric model as a whole or its part can be defined 
by a solid geometry or a FE mesh. In the case of pure geometry, the model is 
discretized by T3D called from MIDAS. However, most of the subsequently listed 
features may be applied to both representations. 

The raw data loaded by MIDAS are parsed in order to build an internal ob- 
ject structure representing the analysed model. On top of that, the complete 
topological connectivity of the model is internally assembled in such a way that 
each geometric element (point, edge, face, cell) is aware of other elements with 
shared vertices. The structured data can be analysed, modified or refined in var- 
ious ways, all done by intrinsic MIDAS's features. These are, for instance, the 
mesh quality control, searching and merging identical nodes and finite elements, 
detection and removal of elements of zero dimensions, localisation and elimina- 
tion of domains separated from the main body, detection of unsupported nodes 
of local kinematic mechanisms, parallel computing support, etc. Multiple inde- 
pendent non-conforming meshes can be connected utilising hanging nodes or rigid 
arms, thus for instance, the effect of reinforcement bars can be integrated in parent 
meshes. Moreover, eccentric joints of beam elements are also allowed through rigid 
arms, where the perpendicular distances between the beam axes are found either 
automatically or fed from the input. 

Structural analysis output data are adjusted to conform with post-processing 
and visualisation. In particular, we plot the cross-section resistance ratio u^i that 
ranges from to oc and has the elastic-plastic threshold Uei^um = 1- Its evaluation 
is based on the von Mises yield criterion 

f{a,k) = ^2-k = 0. (1) 

Assuming the equivalent stress in the form 

(^eq = \/3J2, (2) 

we can write 

Uel = (Jeq/Ry, (3) 

where Ry denotes the yield stress and J2 is the second invariant of the stress de- 
viator [21]. It is obvious that u^i < 1 indicates beams loaded in elastic regime 
and values greater then 1.0 those witch are developing inadmissible plastic zones. 
We are fully aware of this indicator being rather provisional, especially for mate- 
rials of anisotropic strength, however, it provides us with an instant and sufficient 
information on the overall stress distribution in the entire model. 

Due to the license regulations covering the source code, MIDAS can be easily 
extended according to additional needs, e.g. when solving non-standard problems 



with complex geometry and topology. Section 5.3 MIDAS is the principal ingre- 



dient of the proposed methodology, as it integrates all the remaining components 
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Figure 5: Demonstration visual program with DONKEY components. 



together. It is a surrogate for the structural engineer's expertise, thereby allowing 
to reduce his/her personal involvement with a post-processing kit. 

In an ideal situation, architectural and structural models are identical and 
MIDAS converts the data received from the modeller directly to the FE solver. In 
particular, it selects material characteristics from the database, assigns them to 
the finite elements, prescribes required loads, kinematic constraints and produces 
the OOFEM input file. In more complicated situations, the model can be refined 
by making use of any of the MIDAS 's features mentioned above. 

12. DONKEY 

The graphical algorithm editor Grasshopper ^22j, closely integrated with the 
NURBS-based 3D modelling tool Rhinoceros, was chosen as the coding frame- 
work of the plug-in DONKEY [HI [15]. Grasshopper is a visual programming tool 
for procedural modelling popular among academics and professionals. It allows 
designers to generate simple geometries as easily as the awe-inspiring ones, still 
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Figure 6: Cantilever and cross-section resistance ratio drawn on deformed shape. 

preserving possibility of interactive modifications. Programs are created by drag- 
ging components with particular functionality onto a canvas. The outputs of these 
components are then connected to inputs of subsequent components. In this en- 
vironment, DONKEY is accessible as a set of components in a separate tab of 
Grasshopper's menu, see Fig. [5| Properties of DONKEY components highlighted 
in Fig. [5] are demonstrated on the example of a cantilever of 1000 mm in length 
and circular cross-section of 20 mm in diameter, being subject to the vertical force 
of 264 N, corresponding to 27 kg, acting at the unconstrained tip. 

In the first step, a user creates a geometric model, appropriate for structural 
analysis, by using a fully automated tool (algorithmic architecture) or a standard 
drawing procedure (human input based CAD layout). In this particular example, 
a single line and its end points are obtained by Grasshoppers' built-in functions 
Line and End Points, notice the group Geometric model in Fig. [5} Within 
the second step, each of the entities is provided with the information necessary 
for numerical analysis, thereby the solid geometry becoming structural model. In 
particular, circular cross-section and steel material was assigned to the line by 
components Profile and Steel. Next, a constraint (all displacements and ro- 
tational degrees of freedom constrained by default) is applied to one, component 
Support, and a force load to the opposite end point of the beam, component 
Load. The model, component Model is exported to a VTK XML file. Tab. [l| 
and sent to MIDAS, component Analysis. It is further discretized by calling 
T3D and analysed in OOFEM. Finally, the cross-section resistance ratio and me- 
chanical quantities such as strains, stresses and displacements can be visualised by 
corresponding components, see the group RESULT interpretation. The screen- 
shot of Rhinoceros view-port captures the structural model and the cross-section 
resistance ratio drawn on the deformed cantilever. Fig. [6| The highest calculated 
resistance ratio is L175, as visible in Fig. [5| component Max. 
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4^3. Exchange data file format 

The flow of data proceeds between DONKEY and MIDAS by means of flies in 
VTK XML format, Tab. [TJ The geometry is deflned through initial pair of data 
blocks proceeded by the POINTS and CELLS keywords. Structural properties 
assigned to geometric elements are stored in POINT_DATA and CELL_DATA 
sections. The unstructured section AppendedData contains generic information 
of the project's name, material speciflcations, cross-section characteristics etc. To 
speed up the data flow, the ASCII is replaced with the binary format and the 
particular flies are stored in virtual memory instead of the hard drive. 

5. Case studies 

The proposed concept of integrated design is illustrated on three case stud- 
ies. These were carried out in a close collaboration involving architectural studio 
FLOW at Faculty of Architecture in CTU in Prague (FA CTU), see [23l [24], 
CUBESPACE studio and a free-lance artist Federico Diaq3 All the contribu- 
tors are engaged in algorithmic architecture featuring complex forms generated by 
computer algorithms that are driven by human-entered aesthetic and functional 
contexts [25]. Since these structures diflFer from traditional ones, it is diflicult to 
reliably assess their mechanical behaviour without the computer-aided structural 
analyses. However, a detailed simulation of their response would be too pro- 
hibitive, considering rather early phases of the projects. The results represent the 
flnal responses of manually (The Leonardo Bridge, Annelida) and automatically 
optimised structures (CDF). 

The flrst of three case studies, Kurilla's Annelida bridge [26], represents a 
heterogeneous geometry composed of shells and girders, which requires a signiflcant 
reduction to become an acceptable structural model. On the contrary, the self- 
supporting Leonardo's bridge is much less complicated and the architectural wire 
model almost coincides with that for structural analysis. Finally, we demonstrate 
the full power of MIDAS interface on the investigation to very complex sculpture. 
Geometric Death Frequency- 141, by Federico Diaz [27j. 

5.1. Annelida 

Annelida bridge exemplifles a complex task whose computational model has 
to be signiflcantly simplifled before the structural analysis execution. The bridge, 
made up of steel as suggested by Lukas Kurillajj [26j , was truncated for demon- 
strative purposes toa44x8xl2m segment. The frame is composed of straight 



^www.cubespace.eu, www.fediaz.com 

^http://www. studioflorian.com/projekty/63-lukas-kurilla-annelida 
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<VTKFile type="PolyData" version="0. 1" byte_order="LittleEndian"> 
<PolyData> 

<Piece Number0fPoints="2" NumberOfLines="l"> 
<Points> 



Number0fComponents="3" f ormat="ascii"> 



ascii' 
ascii' 



type="Int32' 
type="Int32' 



Name="connectivity"> 1 </DataArray> 
Name="of f sets"> 2 </DataArray> 



ascii" type="Int32" Naine="Boundary_Conditions" NumOf Comp="6"> 



ascii" type="Int32" Name="ID_BOUNDARY_CONDITION"> 



<DataArray type="Float32 
0.0 0.0 0.0 
1000.0 0.0 0.0 
</DataArray> 
</Points> 
<Lines> 

<DataArray format' 
<DataArray format' 
</Lines> 
<PointData> 

<DataArray format' 
111111 

</DataArray> 
<DataArray format' 

1 
</DataArray> 
</PointData> 
<CellData> 

<DataArray format' 
<DataArray format' 
</CellData> 
</Piece> 
</PolyData> 
<AppendedData> 

<Characteristics> 

<COMMENT> <item> example - cantilever </item> </COMMENT> 
<CROSS-SECTIONS Number="2"> 

<item> 1 Rectangle width 0.1 height 0.2 refNode y -2 </item> 
<item> 2 Circle width 20.0 </item> 
</CROSS-SECTIONS> 
<MATERIALS Number="l"> 

<item> 1 IsoLinEl E 210.0e+03 nu 0.20 tAlpha 0.000012 density 7850.0e-09 </item> 
</MATERIALS> 
<BOUNDARY_CONDITIONS Number="l"> 

<item> 1 NodalLoad components 6 -264.777 0.0 0.0 0.0 0.0 0.0 </item> 
</BOUNDARY_CONDITIONS> 
</Characteristics> 
</AppendedData> 
</VTKFile> 



ascii' 
ascii' 



type="Int32' 
type="Int32' 



Name="ID_CROSS-SECTION"> 2 </DataArray> 
Name="ID_MATERIAL"> 1 </DataArray> 



Table 1: VTK XML file generated by Donkey. 

and arc tubes of circular cross-sections that form a repeating geometric pattern 
of distorted rectangles with circular openings, Fig. [7^. The frame vertices are re- 
inforced with a pair of concave steel plates of mutual distance equal to the outer 
diameter of the frame tube being aligned with, Fig. [TJd. The structure is supported 
at two pairs of points, each pair located at a single bridge end. The parametrized 
model was generated automatically by means of a single purpose script. These 
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Figure 7: Annelida, a) complex architectural model, b) detail of joint, c) structural model, d) 
structural model of joint. 

parameters were optimised on the basis of the resulting mechanical response. 

As shown in Fig. [7)3, the architectural model has been created including all the 
details specific to design features and omitting any computational simplifications. 
Even with a high performance computer at hand, it would be barely possible to 
generate a mesh of shell or volume finite elements resolving the model in detail, 
see e.g. [6], as the elements in the tube walls would be much smaller then those 
in reinforcing plates. Such a fine FE discretization would result in an excessive 
computational overhead. Furthermore, technicalities, such as connecting the pairs 
of straight frame bars would be difficult within the "detailed discretization" con- 
cept as well, since these are slightly non-parallel thanks to the distorted topology 
of the entire structure. For these reasons the script was modified to generate a 
simplified architectural model where the frame bars reduce to the centroid axes of 
zero cross-sectional area and only a single mid surface represents the twin corner 
haunches. Fig. [7]1. 

Next steps were identical to the previous cantilever example, with the exception 
of FE mesh generation that was executed in Rhinoceros. Identical coordinates were 
prescribed to all nodes at the contact among beam and shell elements and arising 
multiplicities were merged in MIDAS. 
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Figure 8: Sports hall. 




^/^ 



a) 



b) 



Figure 9: Self-supporting arch a) 3d scheme, b) 2d structural model. Dashed lines indicate 
dosed (with) and open (without) variant. 

5.2. The Leonardo Bridge 

By means of the Leonardo bridge, we would like to demonstrate learn and 
form-finding capabilities of the proposed software. The project of a sports halQ 
for up to 300 spectators was designed by Martin Cisaf , Fig. |8| and analysed by 
Karolma Maskova [28], the undergraduate students at FA CTU and Faculty of 
Civil Engineering of CTU in Prague, respectively. The structure is composed 
of fourteen arch sections inspired by Leonardo da Vinci's self-supporting bridge, 
famous for its ingenious simplicity. Fig. |9^. Besides its structural efficiency, the 
system is known for the self-locking joints, which enable fast erection without 
fasteners and easy disassembly. 

The typical arch is of 35 m in length and 13 m in height. It is assembled of 



^www. studioflorian.com/projekty/184-martin-cisar-mestska-sportovni-hala-v-kutne-hore 
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a) 



b) 




Figure 10: Three distinct models with identical setup and displacement - a) 
closed variant and c) closed variant with mobile support. 



open variant, b) 



timber beams rectangular in cross-section, which must resist loading by the self- 
weight (the segments themselves plus the roofs dead load) and standardised weight 
of snow. The structural model, Fig. [9)3, was generated by an algorithm with pa- 
rameters of the number of segments and lengths, and cross-section dimensions of 
individual beams. At the first instance, three distinct models with identical setup 
of design parameters were compared. Fig. [TOJ Although the open variant. Fig. |9] 
and Fig. 10^, is more advantageous from the application point of viewj we ob- 



serve the significant local displacements arise beams adjacent to applied supports, 
compared to its closed counterpart. Fig. [T0) 3. In particular, the maximum total 
displacements and the cross-section resistance ratio are 44 mm and 0.45 for the 
closed variant in contrast to 154 mm and 1.42 for the open one. Another funda- 
mental distinction in overall structural response brings the removal of horizontal 
constraint in one of the supports. Fig. ^. Besides the different shape of fiexural 
curve. Fig. [TOb and Fig. 10:, it is obvious that the self locking mechanism is fully 



allowed only for the arch with the mobile support due to the negative bending 
moment at the top of the arch in Fig. [TOb. 



^simple. wikipedia.0rg/wiki/File:Da_vinci_bridge.jpg, http://www.rlt.com/20101 
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Figure 11: Examples of various shape variants, undeformed and deformed shapes. 

Resulting from the optimisation process above, the closed-form Leonardo scheme 
with fixed horizontal degrees of freedom was selected for the subsequent intuitive 
form- finding process shown in Fig. 11, The arch shape and cross-section dimen- 



sions were adjusted by making use of a parametric script in order to minimise 
displacements and cross-section resistance ratio of beam elements, yielding to the 
optimal shape. 

5.3. Geometric Death Frequency- Uf.! 

The last example is to demonstrate the MIDAS 's capability in application to a 
geometrically complex artwork with a cellular substructure made up of synthetic 
materials, the Geometric death Frequency-14lP](GDF), installed by Federico Diaz 
in the exterior of MASS MoCA (Massachusetts Museum of Contemporary Art) 
exhibition area in 2010 [27] . 



RealFlow 




Figure 12: GDF - flowchart of individual program components. 

GDF represents the 141-st frame of the fluid flow analysis of a certain amount of 
liquid suddenly entering a closed box. The fluid motion was simulated numerically 
by RealFlow [29] and the particular frame was selected as the starting point for the 
subsequent optimisation process based on the static response. The emerging wave- 
like form was spatially filled up with hollow Acrylonitrile Butadiene Styrene (ABS) 
balls of 47 mm in diameter and 1 mm wall thickness by means of single-purpose tool 



Robo.d [30j, Fig. 12, Nearly 420 thousand of balls has been assembled in a regular 
grid and glued together in contact points, thereby forming the self-supporting 
structure. Fig. [13| The huge amount of basic spherical cells made the manual 
fabrication and quality control management of all contact details unfeasible. Hence 
the entire process has been fully robotized. 



^See www. massmoca.org/event -details. php?id=549 for more details; the fabrication process 
can be found at http://vimeo.com/16019145. 
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Due to GDF's structural complexity, the mechanical response to applied loads 
(dead load, snow weight) was difficult in a fully automatic way by making use of 
the basic MIDAS functionality. Moreover, it was required by the author's team to 
implement additional functions for a decision management based on a priory values 
defined by an expert. The output data were therefore simplified to bi-coloured yes- 



no diagrams (beams with exceeded bearing capacity are in black. Fig. 14). 

To speed-up the numerical analysis, only the compact arch-shaped part of the 
structure comprising of about 250 thousands balls was considered as critical. The 
sphere-shell sponge-like composite was transformed into a beam finite element 
mesh with nodes placed in the sphere centres. Thus, the beam elements represent 
hourglass like rotational surfaces made up of two half-spheres connected at their 



poles by a droplet of glue. Figs [13] and [14} Although such a geometry yields 
a variable stiffness, the beams were considered as prismatic and with averaged 
material characteristics. The bearing capacity of the homogenized beams, normal 
and bending stiffness were obtained experimentally by the load test of several 
cantilever girders consisting of ten axially aligned balls. The measured quantities 
were verified by a detailed FE analysis of a three-dimensional model with the 
balls and glue joints precisely resolved. The material parameters of ABS plastic 
and the glue were provided by the manufacturer. Finally, structural supports 
(displacement constraints) of the arch model were applied to all nodes representing 
contact points among spheres and the horizontal base. 

The transformation process of GDF solid representation into the FE model was 
controlled by Robo.d. Despite the fully automated conversion, the raw mesh was 
further validated by MIDAS interface. First, nodal and element duplicities were 
eliminated. The nodes and elements separated from the central body mass (arising 
from separated drops or splashes of liquid) were identified and excluded from the 
analysis and subsequently from the sculpture itself [27j. Loads and boundary 
conditions have been applied after the adjustment. 

The FE model was exported to VTK file and revised visually in order to find 
any major defects owing to the automatic processing (model connectivity, over- 
all geometrical deviations between FE and the solid model, etc.). Next, it was 
controlled once again by the OOFEM preprocessor routines and analysed. The 
resulting mechanical quantities were post-processed by MIDAS and visualised in 



Paraview, Fig. 14, 



The numerical model contained about 800 thousand degrees of freedom. There- 
fore, the iterative IML [31j solver of the global algebraic system with incomplete 
Cholesky preconditioning was used. This solver, however, exhibits poor conver- 
gence for structures with non-uniform stiffness distribution. In this particular case, 
such an inhomogeneity was attributed to elongated splashes of the liquid (dead 
arms). Thus, the function eliminating the arms of 1 to 2 balls in diameter was 
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Figure 13: Geometric Death Frequency- 141, visualisation. 






'Wi 




Figure 14: Arch-shaped part of GDF - cross-section resistance ratio rendered as bi-coloured 
scheme. Black elements indicate values greater than 1.0. 



further implemented to MIDAS. This led to removal of 1 % of finite elements with 
a negligible effect on overall response while reducing the computational time down 
to fractions of the original time. 



Solving the structure, certain floppy spots had been detected, Fig. 14, The 
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shape evolution then proceeded to choosing yet another frame of the fluid stream 
and either incorporating or removing some ABS ceUs in appropriate regions. This 
was repeated several times until the frame 141 and its optimal shape appeared. 

6. Conclusions 

This article is devoted to the initial component of the integrated design of 
geometrically complex structures, in particular, to the simulation of a mechanical 
response in the conceptual phase of architectural design. It aims at maximum 
possible automation of structural behaviour assessment in the early stages of the 
design and results in economic and reliable exploration of designer's creativity. 
A simple, though effective methodology based on an open source interface that 
allows for interconnecting existing computer aided design and structural analysis 
engineering tools was introduced. Based on three illustrative case studies, it can 
be conjectured that: 

• if the architectural model is created with respect to a subsequent structural 
analysis, the proposed process is robust and free of further user intervention; 

• in the case of a complex model, certain simpliflcations to the model are 
required, however, the structural analysis can be still performed without the 
need for structural engineer's interventions; 

• on the contrary, collaboration with experts in structural analysis, numerical 
methods and programming is necessary when solving extraordinary and/or 
very large structures; 

• signiflcant time savings in communication between structural engineers and 
architects were achieved when solving all three benchmarks, no matter the 
complexity. For example, 20 modiflcations of the Geometric Death Frequency- 
141 model was made within 14 days. 

Finally, let us emphasise that our aim is not to replace a detailed structural as- 
sessment up to the extent required in the advanced stages of the project (building 
certiflcate and/or operating documentation) but to provide architects, designers 
and artists with a simple tool assisting in better understanding of structural be- 
haviour. 
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